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ABSTRACT: A systematic approach is followed in order to identify realistic D-brane models 
at toric del Pezzo singularities. Requiring quark and lepton spectrum and Yukawas from 
D3 branes and massless hypercharge, we are led to Pati-Salam extensions of the Standard 
J> ■ Model. Hierarchies of masses, flavour mixings and control of couplings select higher order del 

Pezzo singularities, minimising the Higgs sector prefers toric del Pezzos with 0IP3 providing 
. the most successful compromise. Then a supersymmetric local string model is presented 

with the following properties at low energies: (i) the MSSM spectrum plus a local B — L 
gauge field or additional Higgs fields depending on the breaking pattern, (ii) a realistic 
hierarchy of quark and lepton masses and (iii) realistic flavour mixing between quark and 
lepton families with computable CKM and PMNS matrices, and CP violation consistent 
with observations. In this construction, kinetic terms are diagonal and under calculational 
^ ■ control suppressing standard FCNC contributions. Proton decay operators of dimension 

4, 5, 6 are suppressed, and gauge couplings can unify depending on the breaking scales from 
string scales at energies in the range 10 12 — 10 16 GeV, consistent with TeV soft-masses 
from moduli mediated supersymmetry breaking. The GUT scale model corresponds to D3 
branes at 0IP3 with two copies of the Pati-Salam gauge symmetry SU(4) x SU (2)^ x SU (2) l . 
D— brane instantons generate a non-vanishing \i— term. Right handed sneutrinos can break 
the B — L symmetry and induce a see-saw mechanism of neutrino masses and R-parity 
violating operators with observable low-energy implications. 
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1. Introduction 

Constructing the Standard Model or a realistic extension within string theory is one of the 
biggest problems in string phenomenology. The challenge is magnified due to two competing 
facts: the large degeneracy of string models and the many experimentally verified properties 
of the Standard Model and evidence beyond. 

On the first issue, a large amount of discussion has been concentrated on the huge 
degeneracy of string vacua in four-dimensions, given by the number of Calabi-Yau com- 
pactifications, the choice of Standard Model embeddings, fluxes, etc. This lack of unique- 
ness complicates the extraction of concrete model independent predictions of string theory 
that can be subject to experimental test (besides the standard 'predictions' concerning the 
existence of gravity and other interactions, dilaton and axion-like fields, moduli, extra spa- 
tial dimensions and supersymmetry but without setting the relevant mass scales). Much 
discussion has been dedicated to this lack of uniqueness leaving sometimes the 'existence' 
of realistic solutions almost as a non- issue. 

This second issue is actually major since contrary to non-stringy model building, a 
successful string model has to be realistic in all aspects, that means it has to be consistent 
with all observational constraints in both high energy physics and cosmology. In particular 

• It has to give rise to the gauge structure of the Standard Model with all the observed 
interactions, the three families of matter fields and at least one Higgs field or an 
alternative mechanism of symmetry breaking. 

• It has to explain the hierarchy of masses of quarks and leptons, including neutrino 
masses. 

• It has to accommodate all flavour issues such as the right mixing angles in the CKM 
and PMNS matrices and right amount of CP-violation, preventing the existence of 
unobserved FCNC. 

• It has to explain the hierarchy of observed gauge couplings either by a consistent 
unification at the appropriate computable scale or otherwise. 

• The proton has to be stable enough to be consistent with observations while a concrete 
mechanism for baryogenesis has to be present. 

• It has to account for the dark matter of the universe. In particular it has to explain 
why extra matter fields do not cause cosmological problems, such as over-closing the 
universe. 

• It has to account for inflation or other alternative early universe mechanism that 
addresses the same problems and gives rise to the density perturbations of the CMB. 

• It has to address the hierarchy problem by low-energy supersymmetry or otherwise 
and has to determine dynamically all the relevant dimensional scales, such as the size 
and shape of all extra dimensions that determine the string and Kaluza-Klein scales 
in terms of the Planck mass, but also the electroweak scale. 
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• It has to have a mechanism to address the dark energy problem. 

The important point that is usually overlooked is that even if a string model successfully 
accounts for most of the observable constraints, it takes just one of these conditions not to 
be satisfied to rule out the model. Due to this problem, much effort has been concentrated 
in identifying mechanisms in classes of models rather than in searching for explicit models. 
But this does not fully spare us from the need to have explicit constructions of realistic 
models in order to be consistent with the standard claim that string theory is the best 
candidate for a unified theory of nature and to at least serve as potential benchmark points 
for physics beyond the standard model. 

In the past decades there has been significant progress on building chiral string models 
of Particle Physics [1]. However to date there is no single compelling model that can 
incorporate the matter of the Standard Model, the spectrum of masses and the hierarchies 
in mixing among generations of quarks and leptons and the other points mentioned above. 
This is mainly due to lack of control over bulk geometries. This problem can be partially 
avoided by focusing on local model building. There, the bulk geometry can be locally 
controlled in the limit of decoupled gravity as for example in models built in intersecting 
D-brane models [2, 3], F-theory [4-9] or D-branes at singularities [10-14], the latter of 
which we focus on in this article. In such models the Standard Model matter content 
and couplings are completely determined by the local geometry and model building can be 
treated independently from other problems such as moduli stabilisation, dark energy, or 
supersymmetry breaking. This approach is known as bottom-up model building, which was 
initiated in [10]. 

Local models are also perfectly suited to be incorporated in the success story of moduli 
stabilisation in type IIB string theory. This has allowed to determine all the relevant mass 
scales, addresses supersymmetry breaking and even the dark energy problem (although still 
subject to debate, it provides the mechanism that makes calculations neglecting this prob- 
lem meaningful). Of particular interest is the so-called Large volume scenario [15], allowing 
for moduli stabilisation using a' and non-perturbative corrections, since it implements a 
concrete realisation of the modular approach to string model building. Typically, shrinking 
4-cycles, the different 'modules', are responsible for the solutions of different physical prob- 
lems such as moduli stabilisation, supersymmetry breaking, cosmological inflation [16,17], 
and the realisation of the Standard Model. But at least one of the moduli is exponentially 
large and does not host the standard model, implying the standard model lives at branes on 
small cycles and then a local origin to the Standard Model. In this context we are interested 
in building supersymmetric models locally from the bottom-up. There are various ways of 
breaking supersymmetry in the local model, some of which construct a sector geometrically 
connected to the singularity as for instance in [18], and others which rely on supersymmetry 
being broken by a geometrically separated sector such as in the Large volume scenario. The 
characteristic soft-masses for local models have been developed in [19-22]. 

In this article we focus on gauge theories arising from D-branes at toric singularities. 
This guarantees the absence of complex structure moduli in the superpotential. This class 
of models offers not only a very rich structure of interacting unified gauge theories with 
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known superpotentials but also leading order control over the Kahler potential. Both the 
superpotential and Kahler potential are controlled by global symmetries, arising from the 
underlying geometric background. With this knowledge of the gauge theory at hand we 
want to demonstrate that these structures are sufficient for realistic Yukawa couplings both 
in quark and lepton sector. By restricting ourselves to model building in perturbative type 
IIB string theory the matter representations in combination with perturbative realisations of 
the Yukawa couplings restrict our model building to non simple-group grand-unified models 
(i.e. we are restricted to models based not on SU(5) or 50(10) which is OK since this is 
not a requirement for coupling unification in these models). Then from the standard GUT 
models, the largest unifying gauge groups we can utilise for model building are the Pati- 
Salam group, 5£/(4) x SU(2) L x SU(2) R [23,24], trinification SU(3) 3 , etc. The Pati-Salam 
model is also the minimal model that realises hypercharge within the non-abelian gauge 
groups and has only one common Yukawa coupling for each family. From the singularities 
perspective we argue that the third del Pezzo singularity (dP^) is the minimal one for which 
the hierarchy of Yukawas and flavour structure is achievable and it is the maximal that is also 
toric. These properties make dP% an optimal candidate to search for realistic models. We 
will present the minimal model that not only includes the Standard Model matter content 
in the Pati-Salam gauge group but also realises its breakdown to the Standard Model 
spectrum. Further motivations for the Pati-Salam model are mentioned in the plethora of 
brane models based on Pati-Salam gauge groups [25-31]. 

A crucial question after fixing the matter content is whether we can achieve the correct 
hierarchy in couplings for quarks and leptons. We obtain the correct quark flavour mixing as 
discussed in [32]. The difference in the lepton Yukawa couplings is achieved via a Majorana 
seesaw mechanism, which can be realised by giving a vev to a right-handed sneutrino. 
We show that the lepton flavour mixing can be different from the quark flavour structure. 
Using the right-handed neutrino masses obtained from giving a vev to the sneutrino, we can 
compute the PMNS neutrino mixing matrix. The precise values of the entries of the matrix 
depend on the details of the RH neutrino masses. However, our model unambiguously 
predicts that the mixing angle #13 is greater than zero, in accord with recent results from 
the T2K experiment [33] . The only toric model that can accommodate all criteria including 
all Standard Model matter is based on the third del Pezzo surface (IP3. 

Gauge coupling unification occurs at the string scale M s or the winding scale RM S 
depending on the appearance of threshold effects [34,35]. The common gauge coupling 
at the unification scale, given by the dilaton, arises from the common geometric origin at 
the singularity for all different gauge factors therefore a simple group GUT is not required 
for unification. Depending on the breaking scale we find various possible unification scales 
ranging from an intermediate string scale at M s ~ 10 12 GeV to M s ~ 10 16 GeV. This flexi- 
bility in the unification scale is very attractive to allow for various string inspired models of 
supersymmetry breaking, partially requiring a string scale below the usual unification scale 
M s ~ 10 16 GeV. Various interesting low-energy phenomenological extensions to standard 
MSSM physics, including additional U(l) symmetries at low energies, additional SU(2) 
vector bosons or additional Higgses, are presented and can be roughly within reach for the 
LHC. 
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Proton decay operators are suppressed by U(1)b-l an d anomalous £7(1) symmetries. 
If they are generated it is at the non-perturbative level after the breaking of U(1)b-l- 
We not currently have a dynamical mechanism to achieve the precise vevs responsible for 
breaking to the Standard Model or to minimise the scalar Higgs potential, but we do check 
for D- and F- flatness. A complete discussion of the latter requires more information about 
the embedding into the compact geometry and is hence tied to moduli stabilisation. 

The rest of this paper is organised as follows. We start with reviewing the gauge theory 
arising from del Pezzo singularities, argue for the special properties that make dP% appealing 
for model building and show how to embed the Pati-Salam gauge group in this setup. We 
then show how non-perturbative superpotential couplings induced from Euclidean 3-branes 
(E3) can give rise to the /i— term. In Section 3 we discuss the breakdown of the Standard 
Model gauge groups, the masses of quarks and leptons and mixings. We also comment 
on the absence of proton decay and analyse gauge coupling unification in this model and 
finally in Section 4 we discuss whether the breaking scales can be combined with TeV— scale 
soft-masses. 

2. A Review of the Ingredients 
2.1 Models at del Pezzo singularities 

There is an infinite number of gauge theories from branes at singularities. The largest class 
of singularities where we control the gauge theory are based on toric singularities. 

The simplest examples of such singularities are the del Pezzo surfaces dPo — dP% where 
one cycle collapses to zero size and which, fortunately, allow for interesting low-energy 
physics. Recall that the n th del Pezzo surface dP n (n = 0, . . . , 8) corresponds to P 2 blown- 
up at up to eight points, the first four dPo — dP% being toric. The gauge theories probing 
these singularities can be obtained using dimer techniques (for a review of the del Pezzo 
examples see for example the discussion in [32]). Up to dP^ no complex structure moduli 
appear in the superpotential [36] , the value of complex structure moduli for higher del Pezzos 
depends on the moduli stabilisation scenario, hence rendering the analysis of couplings more 
difficult and less attractive. In addition dP% is sufficient for our purposes, and allowing us 
to stop there. 

Historically local model building started on C 3 /Z 3 = dPo [10]. Although a realistic 
Standard Model-like matter content can be achieved, the mass eigenvalues are found to be 
of the form (0, M, M), rendering the model un-realistic. To change the couplings one can 
either consider non-trivial B-flux threading the cycles of the singularity [37] , which we shall 
not pursue further in this article, or consider singularities which allow for richer couplings 
per se. In [38] it was realised that models based on dP\ lead to a potentially hierarchical 
mass structure of the form (0,m,M). As discussed in [32], the structure of couplings is 
in-sufficient to achieve the hierarchical mixing angles in the CKM matrix, favouring models 
based on d?2 or higher del Pezzos. As we are interested in models where all Standard Model 
fields arise from D3 — D3 states, we need a mechanism to discriminate between couplings 
for quarks and leptons. As discussed later in Section [3.4[ , this requires an extended field 
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content requiring to consider models based on AP3. 1 Hence from now on we focus on models 
based on the dP^ singularity. 

2.2 Models from the dP% singularity 

The third del Pezzo surface dP% features six U(N) gauge groups and there are four toric 
phases connected via Seiberg duality [39]. Here we choose the toric phase that is connected 
simply via Higgsing to the lower del Pezzo surfaces, its matter content is shown in Figure [l]. 
Recall that every node corresponds to an U(N) gauge group and an arrow from node A to 
B in the quiver corresponds to a field Xab transforming in the bi-fundamental (Na,Nb) 
under U(Na) X U(Nb)- The superpotential is given by 



W, 



dP 3 



-X12Y31Z23 — X^Y^Z^q + X^YsiZiA^p- + X12Y25Z56- 6 



A 



1 Y V 7 * 42 Y V 7 P53$61^42 
+-^36^64^23-^ ^36-125^14 ^3 
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42 

A 



A ^4 
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-X\2 



^12 — 





where A is an appropriate UV cutoff. 2 



A 



(2.1) 




Figure 1: The dP 3 quiver with fields labelled in correspondence to the superpotential in Equation E.ll. 



The structure of the superpotential is entirely fixed by the geometry, which - from a field 
theoretical point of view - can be seen as a global E3 = SU(2) x SU(3) symmetry and 
an additional R-symmetry [41]. The charges of the fields under the global symmetries are 



x lt would be very interesting to see whether this constraint can rule out models based on lower del 
Pezzos whose coupling structure is changed with B-fluxes. Since this option 'only' allows for changing the 
coefficient in front of every coupling, this option does only change the Pati-Salam Yukawa coupling but not 
Standard Model Yukawa couplings for quarks and leptons separately. 

2 Following [40] we assume that within the low energy A/" = 1 supergravity, A is Mp due to holomorphy, 
and the actual physical suppression scale is determined by terms in the Kahler potential. 
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summarised in Table |j. The superpotential can be determined from invariance under these 
symmetries. 



Fields 


SU(2) x 5i7(3) 


U(1)r 


0^36 > ^25> ^14, ^42, 053, $6l) 
PQs' ^56) ^64) 
( Z 23: *31> X ll) 


(2,3) 
(1,3) 
(1,3) 


1/3 
2/3 
2/3 



Table 1: Charges of D3 — D3 states under global symmetries, taken from [41]. The last two rows are 
written separately, since gauge invariance forbids a coupling between them. 

Choosing two copies of the Pati-Salam gauge group as shown in Figure |2[ we have that 
all Standard Model fields are given by D3 — D3 states. 




Fi gure 2: The quiver diagram of the Pati-Salam realisation on dPz. The fields ($61,^42,^53) will be 
responsible for gauge group breaking, combining the two copies of every gauge group. The other fields are 
Standard Model fields including three generation of Higgs fields. 



In particular we identify the fields Q l L = (X45, Y25, Z23) as left-handed quarks (and leptons), 
Q l R = (-Xgg, Y31, Zj|) as right-handed quarks (and leptons) and the matrix as the Yukawa 
matrix. The fields Hi = (Xi2,l64, ^14) are Higgs fields. With this identification we can 
re-write the superpotential in the usual form 



WD3D3 



-H1Q2Q3 ~ Q1H2Q3 + QiQ^Hs^- + HiQ 2 Q R —^- 



+QiH 2 Q3—^- - Q1Q2H3 



P53 

A 

P53$61^42 



A 3 




/ 





Hz p -f -H 2 \ 








^3 A 3 


H^ 






V 


H 2 - r 


-Hi J 







(2.2) 



which we shall use from now on. 
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2.3 Kahler potential 

In a purely local model the global symmetries or isometries restrict the Kahler potential to 
a diagonal form. However this is no longer guaranteed in an effective supergravity setup 
where these global symmetries are broken by volume dependent effects. The precise scale 
at which they are broken is not known since control over the Kahler potential is very 
limited. However, the extended gauge structure of (IP3 provides us with new control over 
the appearance of the matter fields in the Kahler potential. 3 

Since there are only single fields between any two gauge groups, the constraint of gauge 
invariance allows only flavour diagonal terms in the Kahler potential at leading order. Since 
all Standard Model matter fields have the same geometric origin, they will have a common 
overall factor depending on the Kahler moduli. To leading order in the 1/V large volume 
expansion, we can write 

^matter D ""^/a'^ {Ql,rQl,R + + $61^61 + ^42*42 + P53P5 3 ) , (2.3) 

where /(t s ,t&) is a function of the small modulus being suppressed by higher inverse pow- 
ers in the volume. 4 This structure also ensures that the soft-masses induced via moduli 
mediation in a large volume setup are flavour-diagonal, satisfying the conditions presented 
in [46]. In particular, problems for flavour changing neutral currents arising from a scalar 
mass matrix that is not proportional to the Kahler metric K a b are absent as long as the 
gauge symmetries are unbroken [45]. We also note that terms like H3H3 + H3H3, utilised 
in the Giudice-Masiero mechanism, are forbidden by the anomalous U(l) symmetries. 

2.4 Anomaly cancellation and D7 branes 

In order to cancel anomalies for the given choice of D3 brane gauge groups, it is necessary 
to have D7 branes that lead to additional D7 — D3 states. The spectrum of D7 — D3 
states for general toric singularities was developed in [47], leading to one pair of D3 — D7 
and D7 — D3 states for every D3 — D3 state. For every 33 state $3^3^ , there exists a 
supersymmetric 7-brane giving a (7i) fundamental and a (7j) anti-fundamental with the 
Yukawa coupling $3^ ,(7i)(7j). The most general (IP3 quiver including D7 branes is shown 
in Figure ||. As discussed for example in [48], the cancellation of anomalies is checked by 
counting the arrows coming in/out to any node weighted by the rank of the gauge group 
they originate/end. In our (Pati-Salam) 2 model we find the following constraints for the 
D7 rank gauge groups: 

77112 = 4 - m 4 + m 5 + m 6 , 
mil = 2 — mi + ^2 + m A j 
rriio = 4 - mi + m 3 + m 7 , 
m-9 = tti-2 — m3 + m5 + mj — m-j , 

m$ = -2 + m 2 - m 3 + m 5 . (2.4) 

3 Similar use of multiple (7(1) symmetries has also been made in recent heterotic model building [42]. 

4 The precise structure of this modulus weight is not known but can be estimated in various setups [43,44]. 
The structure of the moduli weights is of particular importance for the soft masses as discussed in [20,45], 
and in the context of guaranteeing F-flatness. 
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Solutions with mi > are physically relevant. Note that although we have in principle 
six anomaly cancellation conditions, one of them is trivially satisfied when the other five 
are satisfied. In section |3.2| (cf. Figure we specify a choice for the D7 gauge groups 
compatible with the breakdown to the Standard Model gauge symmetries. We note that 
the additional D7 — D3 states can decouple from the low-energy effective action, since 
D7 — D7 interactions in the bulk away from the local construction can give rise to large 
masses. 



mil 




m8 



Fi gure 3: The dPz quiver with D7 branes included. Indicated in red are the fields that we vev to break 
the (7(4) x (7(4) symmetry, the D3 — D7 states need vevving to satisfy D-term equation. Similarly the states 
which break the left and right U (2) symmetries are indicated in purple and orange. The states highlighted 
in green are D7 states that we could vev in accordance with the D-term equations to give large masses to 
some of the Higgs fields. 

2.5 Anomalous U(l)s 

The U(l) factors contained in U(N) factors appear differently in the low-energy description 
depending on their origin. We distinguish between: 

1. U{1) symmetries within non-abelian symmetries such as the diagonal generators in 
SU(N). These U(l) symmetries remain massless at low-energies. 
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2. Anomalous £7(1) symmetries get string scale masses. 

3. Non-anomalous £7(1) symmetries remain massless in the non-compact model. How- 
ever, once embedded into a compact model they remain only massless if the associated 
cycle becomes trivial in the bulk, [49]. Otherwise they receive masses via the Green- 
Schwarz mechanism. 

As shown in [38] the masses for non- anomalous £7(1 )s upon compactification are suppressed 
compared to the masses for anomalous £7(1) symmetries and are given by the mass for KK- 
modes Mkk ~ 1/V 2 / 3 . In the philosophy of bottom- up model building we demand that 
all cycles associated to non-anomalous £7(1) symmetries are non-trivial in the bulk. We 
then are left with only non-abelian SU(N) symmetries. The now massive £7(1) symmetries 
remain as global symmetries in the low-energy spectrum. 

2.6 E3 branes and non-perturbative effects 

Non-perturbative contributions to the superpotential can be induced from Euclidean E3 
branes wrapped on 4-cycles passing through the singularities in analogy to D7 branes. 
Instanton induced superpotentials for branes at singularities with E3 branes have been 
studied in [50] (for a review of the subject see [51]). The reader is referred to this reference 
for a derivation of the conditions required for a sufficient zero mode structure. We will 
utilise these effects to generate fi— terms for the Higgs fields. 

For our purpose the additional zero modes between E3 branes and D3/D7 branes are 
of interest. The spectrum and couplings of these zero modes are the same as if there were 
a D7 brane wrapping the E3-wrapped cycle. 

The presence of the desired couplings however does not imply that a non-perturbative 
contribution to the superpotential is generated, one has to ensure the correct number of 
fermion zero modes. The classical action for the instanton is given by a global piece de- 
pending on the 4-cycle volume they are wrapping, a local piece depending on the twisted 
moduli at the singularity and a piece coming from couplings to bi-fundamental fields 



where rj x denotes the corresponding zero mode with the E3 and D3 brane, (j)^ are twisted 
closed string moduli which we can neglect from now on, <3? r denotes a 33 state, and is 
a coupling constant, for our purposes mainly indicating whether the coupling exists or not. 
In addition we have couplings among D3 — D7, E3 — D3, and D7 — E3 zero modes of the 
general form 



There can be further zero modes involving D7 — D7 states depending on the boundary 
conditions at infinity which we shall not need in the further discussion. We need to in- 
tegrate over the fermionic zero modes which gives rise to the following non-perturbative 
superpotential 




(2.5) 



VE3D7&37VE3D3 ■ 



(2.6) 




(2.7) 
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where fp refers to the E3-D7 zero modes. In order to give a non-vanishing contribution 
every zero mode has to appear precisely once. The scenario presented below for the //— term 
should be seen as to show how such a contribution can arise, keeping in mind that depending 
on the boundary conditions at infinity there could be more zero modes, which can change 
the contribution by additional vevs appearing in the //— term. 

2.6.1 The //-term 

For a viable low-energy Higgs sector we need to have a fi— term for the Higgs fields fiH u .Hd- 
In a model with left-right extension the Higgs field H including both H u and Hd fields, the 
//—term is conveniently rewritten as fiH.H where the product H.H is to be understood as 
a contraction of the ST/ (2) indices with e a p. 

In models from branes at singularities the Higgs fields are additionally charged under 
anomalous U(l) symmetries which forbid this coupling perturbatively. We remind the 
reader here about the arguments presented in [50] on how to generate this coupling non- 
perturbatively via stringy instantons from E3-branes. To project down to the required 
number of two fermionic zero modes we have to require that either the E3 brane is mapped 
to itself by an orientifold projection away from the singularity and that the singularity is 
mapped to a mirror image of itself or alternatively that a single D7 brane is wrapped on 
that cycle [52]. In the bottom- up philosophy we are not looking for a concrete geometric 
realisation of the compact setup, but we require that anomalies are cancelled locally. 

The Higgs field H transforms as (2r,2~l) and there are the following charged zero modes 
in this scenario between the D3 and E3 brane: 

V L = (1,2 L ,1) and rj R = (2 R , 1, 1) , (2.8) 

where the last charge denotes the charge under the E3 brane gauge group. 
In this setup the following couplings are induced 

Charged = J <H V R . (2.9) 
Now we can integrate over the charged zero modes as follows 

Schargcd = j ' dq^dr^^Hij-nf (2.10) 
= e ac e bd H ab H cd . (2.11) 

This then leads to the following non-perturbative contribution to the superpotential as 
outlined above 

W np = Ae~ aTs H.H, (2.12) 

where T s denotes the chiral superfield associated to the 4-cycle volume which the E3-brane 
is wrapping and A denotes a constant depending on complex structure moduli. Note that 
this allows us to obtain different vevs for the Higgs fields due to different cycles the E3- 
branes are wrapping. Note that we do not have any mixing among Higgs field generations 
through the fj, term in the superpotential. 
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2.7 Neutrino see-saw mechanism 



At tree- level there is no Majorana mass term due to the additional U(l) symmetries of the 
singularity. There are two mechanisms which could generate Majorana neutrino masses: 

1. A vev for the right-handed sneutrino can give rise to Majorana neutrino masses as 



for instance in [53-59]. We postpone the discussion of this mechanism to section 3.3 
since it depends on the structure of vevs we demand for a viable flavour structure. 

2. D-brane instantons could lead to the desired Majorana masses. 

In Appendix ^| we show that the second possibility cannot be realised in the model 
presented in this paper. The neutrino mass term is forbidden by the U(1)b—l gauge sym- 
metry, which will not be broken by non-perturbative effects. However, if we were to break 
this symmetry by giving a vev to a 37 string state, the symmetry would become an effec- 
tive global symmetry which could be broken by the string instantons. After doing this, 
both the U(1)b-L and U{1)y groups mix with the abelian U(l) factors from the D7 gauge 
groups. This in turn implies that the Kahler modulus is charged under these symmetries, 
and generates a Stiickelberg mass for both groups. There is therefore a tension between 
having a non-perturbatively generated mass for the right-handed neutrinos and keeping 
hypercharge unbroken. By extending the D7 sector we have found it possible to generate 
a mass term for the neutrinos; however such a term comes from a higher dimensional op- 
erator leading to un-physically low neutrino masses. Of course, it still may be possible to 
have viable instantonic neutrino masses in models based on different gauge groups, or at 
different, possibly non-toric singularities. The construction of models at toric singularities 
with non-perturbative neutrino masses is an interesting challenge. 



3. The Standard Model in dP 3 

After discussing all the ingredients of our model, let us now focus on how the Standard 
Model matter content and couplings can be accommodated in the third del Pezzo surface 
dP%. We start by discussing how the gauge structure in dP 3 can be broken to the Standard 
Model gauge group. We specify the D7 sector and determine the matter content below the 
breaking scale. After that we discuss the phenomenology in terms of the masses for quarks 
and leptons, the flavour mixing, proton decay, the number of parameters we tune and gauge 
coupling unification. 

3.1 Breakdown to the Standard Model 

The breaking of gauge symmetries can be achieved by resolving the singularity slightly, 
which corresponds in the gauge theory to non- vanishing Fl-terms. To satisfy the D-term 
equations then requires certain fields to acquire vevs. Here we only check that our vev 
assignments are actually flat directions of the potential without completely minimising the 
potential. This is due to the lack of knowledge of how to stabilise the moduli associated to 
the slight resolution of the singularity. In a complete model, the 4-cycle volume associated 
to the Standard Model singularity corresponds to a modulus and needs to be stabilised 



- 12 - 



as well. In local D-brane models at del-Pezzo singularities the singular point corresponds 
to a zero 4-cycle volume. For the del Pezzo 4-cycle the usual non-perturbative effects 
cannot be used for stabilisation [60], the leading contribution then arises from D-term 
potentials associated to anomalous £7(1) symmetries which the 4-cycle volume is charged 
under. Having plenty of anomalous U (1) symmetries associated with the singularity it is 
natural to assume that the del Pezzo 4-cycle modulus will be charged under them. At 
tree-level, this implies a stabilisation in the singular regime, explicitly driving the 4-cycle 
volume to zero size. Higher order effects then will become important to determine the size 
of this breaking. 

A detailed analysis of the stabilisation is beyond the scope of this paper and we only 
pursue the search for interesting flat directions at this stage. 

From now we allow ourselves to set the scale of breaking down to the MSSM to be any 
scale below the string scale. Whether we can break directly to the Standard Model or in 
an intermediate step to a left-right model is dealt with in section |3.7[ 

To break £{7(4) we use a method that appeared in the context of SU(5) GUT models 
[61], starting with two copies of £{7(4) and breaking it to diagonal gauge groups with 
bi- fundamental matter transforming as (4,4). Our approach differs from [26] through the 
presence of only one bi-fundamental field between the SU (4) factors, possible through the 
appropriate inclusion of D3 — D7 states. In the dP^ model, we have precisely one field 
transforming as (4,4). Giving P53 a vev allows for the following breakdown: 



{p. 



53/2 



/m \ 

vi 

vi 

\ v 2 ) 



■ {7(4) a x U{A) b -> £[/(3) diag x {7(1) 4 X U(l) d , (3.1) 



where {7(1)4 refers to the {7(1) in £{7(4) and U(l) d is the diagonal {7(1) of the two {7(1) 
factors from {7(4). The D-flat conditions for this vev can be guaranteed by the presence 
of D3 — D7 states with appropriate vev as already mentioned in Equation C.2. The two 
U (2) l gauge groups are broken to the diagonal subgroup by vevving ^42 as follows: 



(^42) 




(3.2) 



This vev satisfies the SU{2) D-term equations by itself. The ?7(1) equations can be satisfied 
by the appropriate Fl-term. The two f7(2)jj factors can either be broken to the diagonal 
gauge group or directly down to {7(1) x {7(1) via the following vev 



61 




(3.3) 



As in the case for the breaking of the U(4) factors, the D-flatness can be guaranteed by the 
presence of D3 — D7 states as in Equation |C.3| . Further details on the breaking including 
a discussion of the masses associated to these breaking fields and how sufficient F-flatness 
is achieved can be found in Appendix |A[ 
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3.2 The model at low-energies 



The presence of appropriate D-instanton effects as discussed in previous sections fixes our 
choice of D7 gauge groups almost entirely. The only flexible choices are rriQ and < mg < 2. 
The most general choice of D7 gauge groups in this context is shown in figure f|. We choose 
the ranks to satisfy the anomaly cancellation condition ( |2.4D and still to allow for consistent 
breakdown of the Pati-Salam gauge group to the Standard Model. 5 




Figure 4: The D3 — D3 and D3 — D7 spectrum in the Pati-Salam model. The ranks for the D7 branes 
are chosen to satisfy the anomaly constraints and allow for breaking to the Standard Model gauge group. 



As discussed in section ^Jj, we assume that all abelian U(l) factors of the D3 gauge groups 
become massive. We identify the non-anomalous linear combination of Tjfy/gs and [7(1)4, 
the diagonal U(l) in SU(4), as hypercharge. 

3 and |3.3| lead to the D3 — D3 spectrum after symmetry 



The vevs in equations [3T 
breaking shown in Table 0. 



5 This restriction is not possible when considering neutrino Majorana masses induced via a vev for the 

A model without D7 branes needed for anomaly cancellation is 



sneutrino as discussed in section 3.3 
presented in Appendix |e| 
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Table 2: The D2> — D3 spectrum after symmetry breaking. For correct hypercharge assignments we need 
in our conventions a = —1/6 and k = 1/2. 



U(l) x is generated by T 15 



5(7(4) 



-,SU(2) r 



. This additional f(l) symmetry, which is just the 



difference between B — L and hypercharge, can be broken dynamically below the string scale 
for instance utilising the right-handed sneutrino as for instance in [53-59], which we discuss 



in section p,3[ Whether such a breaking occurs would be subject to a detailed running of 
the soft-masses from the high-scale which is beyond the scope of this article. 

3.3 Majorana masses via sneutrino vev and R-parity violation 

The alternative to non-perturbatively generated Majorana neutrino masses is to give a vev 
to the right handed sneutrino which was discussed for example by [53-58]. In principle it is 
possible to generate the vev for the sneutrino radiatively [53-56] using the running of soft 
supersymmetric parameters and hence avoiding the problem of D-flatness. As discussed in 
this literature the sneutrino vev generates a mass for the right-handed neutrino of order 
the breaking scale and breaks U(1)b-l x ^(1) to hypercharge, breaking the additional 
U(l) at low energies and leaving only the Standard Model gauge groups. As such this 
option seems very desirable, however it can lead to dangerous R-parity violating operators 
depending on the breaking scale and the right-handed neutrino masses. As detailed above 
the latter are essentially given by the up-quark Yukawa couplings, fixing the breaking scale 
to the range 10 11 - 10 13 GeV. Let us illustrate how we can achieve the absence of R-parity 
violating operators. As discussed in the previous section, we are interested in a solution 
where Hf has a large supersymmetric mass. Looking at (P 2 ) / Owe then have the following 
potentially dangerous coupling 

V2 



M3#3 



A 



(3.4) 



where the first term is a standard fi term for the H3 Higgs generation and the last term is 
a bi-linear R-parity violating term oc H U L. This term can be rotated away by the following 
field re-definition 



m 



m 



V2W2) 



L\ and L l = L\ + 



^2(^2) 



m 



In terms of the rotated superfields and dropping all primes we obtain 



W = Wmssm 



A A 3 



-L 2 Liei 



1)2(1*2) fi# L 

V2 M«l • 



A A 3 



(3.5) 



(3.6) 
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Thus, we have generated lepton number violating terms but not any baryon number gen- 
erating ones. Constraints from baryogenesis [62] restrict both couplings to be smaller than 
1(T 7 . The absolute size of these vevs will not be restricted by requiring realistic hierarchi- 
cal flavour mixings, we still can tune ui^/A, only their ratio will be restricted (cf. Equa- 
tion 3.19). With this freedom at hand, we can clearly satisfy the bound R-parity violating 
operators. Note also that we are attempting to restrict the smallest non- vanishing Yukawa 
coupling which is already heavily suppressed compared to the other Yukawa couplings. 



3.4 Masses and flavour mixing 

The superpotential after breaking becomes: 



w 



+ 



+ 





TTU4> 

U 2 A 



TTU VI 
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fl 3T 



Tjd± 

U 2 A 



TTU V2<t>"4> 

n 2 A 



TTd± 
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TTU V2 

U 3 A 



udV2 
U 3~A 



+ WD3D7 • 



(3.7) 



The D3 — D7 states associated to the breaking fields /O53, ^42, ^61 get vevs to satisfy D- 
term equations. As previously discussed we can generate Majorana masses for right handed 
neutrinos of the following type 



A\v\v\ + A2V2V2 + A 3 u 3 iy 3 



(3i 



where the Ai are some coefficients which absorb all the sneutrino vevs and other couplings. 
For now appropriate configurations for the /i— term are assumed to be present but are not 
detailed here to focus the analysis on the flavour physics in the quark and lepton sector 
whose structure we now would like to analyse. 



Quark sector 

The phenomenology of the quark sector is essentially described by the left-right model 
in [32]. To get a difference between the down and up- type Yukawa couplings we give 
the Hf field a large supersymmetric mass via appropriate D7 vevs. To get the correct 
hierarchical suppression in the CKM matrix we assume the following scaling between the 
Higgs vevs 



Hf 



A 



e and 



(3.9) 
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All contributions, including vl/42/A, are sub-leading. Then the CKM matrix is approxi- 
mately 

/ 1 e e 3 \ 



|Vckm| 



e 1 e 2 
e 3 e 2 1 



(3.10) 



V e e 1 / 

The expression for the mass eigenvalues was presented in [32] and in the above approxima- 
tion they lead to the following hierarchies of masses: 

«) 2 = (0, |fljf + |^| 2 %£, Itffj 2 + l^l 2 ^) , (3.11) 

«) 2 = (0, |# 2 T%^ |tf 2 T + l^T^gT) • (3.12) 

The zero mass for the lightest generation is characteristic of models on toric singulari- 
ties [32]. It can be argued to arise from a global symmetry of the low-energy theory which 
is broken by radiative corrections using bulk effects, leading potentially to realistic masses 
for the lightest generation [63]. 

Lepton sector 

Without a mass term for the right-handed neutrinos we would observe the same Yukawa 
matrices for quarks and leptons at leading order. The difference in the Yukawa matrices 
coming from the breaking of the £7(4) factors is only important at sub-leading order. Re- 
garding the leptons this feature is interesting since the individual down-type quarks and 
lepton masses are not hierarchically different, but it is phenomenologically unrealistic re- 
garding the neutrino and up-type quark masses. The difference might be due to radiative 
corrections. To leading order the lepton masses are given by 

™2 _ f n I ud\2 l^42| 2 i rrrf|2 , | ud |2 \p53 P ^ /o i o\ 

m e L ~ I °> \ H 2\ —£2~i \ H 2\ +\ H ?,\—£r)- (3.13) 

We now turn to the discussion of the seesaw neutrino mass scenario from sneutrino 
vevs. Having Majorana masses for right handed neutrino masses we start with the following 
general mass matrix 

(314) 

This matrix can be brought into block-diagonal form and we then can diagonalise the 
remaining blocks to obtain the eigenvalues of M n .M.\ analytically. We find to leading order 
the following masses for left and right handed neutrinos: 



,2 



2|#f| 2 2|tff| 2 |^42| 2 V (2\H% 2|^| 2 |^ 3 | 2 \ 2 8\H?\ 2 \H%\ 2 \pl 3 \ 2 



A 2 \ \A\\A? J 'V|A 3 | |^ 2 |A 2 ; |.4 2 |A 2 



ml R = (^l 2 ,^! 2 ,^! 2 ). (3.15) 
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We have the seesaw suppression of left-handed neutrino masses as phenomenologically de- 
sired. We note that the left-handed neutrino masses are approximately given by m UL ~ 
m\ ct /m UR . Cosmological constraints imply that m tot = J2i m ^i — 0-5 e V. This implies 
that the right-handed neutrino masses cannot be larger than 0(1O 13 ) GeV. Further con- 
straints from the SuperK experiment imply that at least one of the left-handed neutrinos 
should have a mass of at 0.04 eV. The window for right-handed neutrino masses is then 
approximately 10 11 < m VR < 10 13 GeV. 

Next in line is to estimate the flavour mixing in the lepton sector. The flavour mixing 
is to date known to be present but not as constrained as in the quark sector. The so-called 
PMMS matrix parametrised by three mixing angles can be written as 



C12C13 C13S12 e S13 



(3.16) 



^13 \ 

VpMNS = I -C23S12 - Ci 2 e td Si 3 S 2 3 C12C23 - e td Si 2 Si 3 S 2 3 C13S23 

-Cl2C23e wS Si3 + S12S23 -C23e* 5 Si 2 Sl3 ~ Ci 2 S 2 3 C13C23 J 

Experimentally the mixing angles are constrained as follows [64] 

sin 2 (20i 2 ) = 0.87 ± 0.03, sin 2 (29 23 ) > 0.92, sin 2 (29 13 ) < 0.19 CL = 90% . (3.17) 

In particular with the undetermined phase 5, most entries are not heavily constrained from 
a string model building perspective apart from the fact that there is a large mixing between 
the first and second generation, which we would like to achieve. As in the quark sector 
the mixing between the first and second generation is determined (to leading order) by the 
ratio 

(3.18) 



A 2 Hl 

So far the different vevs of ^53 for lepton and quark sector have not played a role at leading 
order (e.g. masses). Here however, we see that the vevs have to differ in order to predict 
different mixing angles for quarks and leptons. In particular we demand that the ratio of 
both vevs is similar to 

^~e 2 . (3.19) 

The size of this ratio V\jv2 also implies that the mixing angle between the first and 
third generation for leptons is smaller than the corresponding angle for the quarks. This 
can be seen from the scaling of the (1,3) entry in the PMNS and which is given by 

sin 9 l3 PMNS = , , (3.20) 



1 + 



H d n' 2 ' 3> d 
A 2 H$ 



(2) 

which includes a factor of v 2 = p\ 3 in the denominator. The corresponding angle in the 
CKM matrix instead has v\ = and since and since v\ < V2 the mixing angle is smaller for 
leptons. Our model thus does not display tribimaximal mixing. Recent results from the T2K 
experiment [33] have indicated a non-zero value for this angle, 0.08 < sin6»i 3 PMAr5 < 0.27 
being the 90% confidence limits. Since sin^^^ w 0.0035 there is a slight tension between 
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our prediction and the T2K result. However, we expect that due to radiative corrections 
necessary to generate the first family fermion masses that this may also change. Also the 
remaining mixing angle depends on the size of the right-handed neutrino masses which is 
not controlled within the local construction. 

3.5 Counting parameters 

Neglecting the complex phases for now, where we expect some relative alignment arising 
from the scalar potential, we have 9 parameters (6 masses + 3 angles) in the quark sectors 
and 12 parameters (9 masses + 3 angles) in the lepton sector, totalling 21 parameters. In 
our model we choose 15 non-vanishing vevs as summarised in Table |3[ 
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Table 3: A list of parameters that we have to adjust to be non- vanishing in our model to achieve for 
Standard Model masses and mixing angles. 

Our model gives expressions for the 21 SM parameters in terms of these 15 parameters. 
The fact that there are more SM parameters than model parameters implies that there 
must be relations between the expressions for the SM parameters. Some of these are clear; 
for instance m u = = m e = m Ul = 0. There is also a relation between the quark and 
neutrino masses, m UL ~ m u,c,t/ m VR- Finally there are relations between some of the light 
quarks and leptons, which we do not consider robust predictions of our model as they will 
be changed during compactification [63]. 

3.6 Proton decay 

Here we discuss proton decay operators up to dimension six and find that the additional 
(anomalous) U(l) symmetries severely restrict the possible operators. Our analysis is based 
on the list of proton decay operators in Table 3 of [65]. 

Before forbidding potential operators with the anomalous U(l) symmetries, let us recall 
how we can break these U(l) symmetries for which we have two possibilities already used 
in this article: 

1. Using non-perturbative corrections arising for example from E3 branes wrapping a 4- 
cycle intersecting with the singularity (cf. the fj,— term ). In this case a 4-cycle Kahler 
modulus is charged under the U (1) symmetry and can generate, given sufficient zero 
mode structure, the perturbatively forbidden coupling in the superpotential. 

2. A vev for the scalar component of a bi- fundamental field charged under two anomalous 
U(l) symmetries (cf. breaking to the Standard Model gauge groups and a vev for 
the right-handed sneutrino). Couplings breaking these J7(l) symmetries are couplings 
which involve the field breaking the symmetry. No further couplings that can either 
be constructed without the breaking field, or be constructed by integrating out the 
breaking field are induced. Hence for example the dangerous coupling uddv^/K. is 
forbidden by the U(l) symmetries that are broken by the vev of vr and hence also 
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cannot be induced if the right-handed sneutrino obtains a vev. On the contrary a 
higher dimensional Yukawa coupling of the type ^j^QlH u ur remains present in the 
low-energy theory. 

The U(1)r symmetry in the superpotential Eq. [I] restricts the allowed couplings such 
that higher powers of a given superpotential term are forbidden at tree level. The non- 
renormalisation of the superpotential at the perturbative level then guarantees that these 
operators are not induced perturbatively. Hence dangerous operators that one might imag- 
ine involving D3 — D7 or D7 — D7 states are simply absent in the superpotential at the 
perturbative level. On the non-perturbative level and in the Kahler potential, as highlighted 
m Table |, the remaining operators are mostly forbidden by U(1)b-l- The only surviving 
operators are QlQlQlQl an d QrQrQrQr as potential non-perturbatively generated op- 
erators in the superpotential. All other operators in the Kahler potential are forbidden. 
However the operators QlQlQlQl and QrQrQrQr are not generated with the Euclidean 
3-brane setup we choose in our model. In the absence of a global completion of our local 
model, one cannot say whether such operators are induced by other more complicated non- 
perturbative effects. Such effects are beyond the scope of local model building. In summary, 
the R-symmetry and the U(1)b~l forbid proton decay arising from both D- and F-term 
operators sufficiently. 

3.7 Gauge coupling unification 

We now turn to investigating gauge coupling unification in our model. At the high scale 
the gauge coupling for the non-abelian gauge groups is given by the tree-level value of 
the dilaton. With gauge coupling unification implemented at the high scale the question 
to ask is whether we can reproduce the observed values for the inverse gauge coupling 
at the weak scale (o^ 1 = 98.33 ± 0.03, aj 1 = 29.57 ± 0.03, q 3 1 = 8.45 ± 0.05) at M z = 
91.1876±0.0021 GeV [64,66]. In our analysis we restrict ourselves to one loop beta functions, 
neglect threshold effects and restrict ourselves to completely supersymmetric running for 
simplicity. 

The one-loop beta function for a supersymmetric theory with a product gauge group 
G\ x G2 is given by [67] 

P gi = T(R 1 )d(R 2 ) - 3C7 2 (Gi) , (3.21) 

where for an irreducible matrix representation R a , T(R) is defined by Tr[R a R b ] = T(R)5 ab 
(and from now on T(R) = 1/2), d(R) is the dimension of the representation (the dimension 
of the fundamental of SU(N) is d(djv) = A) and 02(G) is the quadratic Casimir of the 
adjoint representation (note that C2(SU(N)) = A). A sum over all field content in the 
above formula is implicit. We use Eq. |3.21| to compute the beta functions for any gauge 
groups that can appear in the breakdown from the dP$ gauge groups at the string scale 
down to the MSSM gauge groups at low energies. Here we keep the most general breaking 
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Table 4: Summary of proton decay operators of dimension 5 and 6 taken from Table 3 of [65]. A denotes 
any possible MSSM field, e represents either the right-handed electron or neutrino in our model. Besides 
being forbidden by U(1)b-l, which for us seems the strongest constraint, anomalous U(l) symmetries 
restrict all couplings. 



pattern possible 

SU(4) a x SU(4) b M A d SU(4) ^ SU(3) x U(1) B -l , (3.22) 

SU(2) La x SU(2) Lb M A d SU(2) L , (3.23) 

517(2)^, x SU(2) Rb M ^ d SU(2) R 17(1) , (3.24) 

where the mass scale Mx denotes the breaking scale of that gauge group. In addition to 
the above breaking scales we allow for two further variable scales, the scale the additional 
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Higgs fields become massive -Mhj ggs and the scale the D7 gauge content becomes massive 
M D7 . 

Recall that the running of the inverse gauge coupling is given by 



a x \M) = a- x \M x ) + ^\og(^- 



(3.25) 



At the breaking scale M x there are matching conditions for the gauge couplings which 
depend on the normalisation of the unbroken group with respect to the larger gauge group 
before the breaking. A detailed analysis of these matching conditions and the running of 
the inverse gauge couplings is given in Appendix [D]. We summarise in Figure || the gauge 
couplings valid in between individual breaking scales. 
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Figure 5: On the left side the breakdown from the dPs gauge groups to U(1)y is shown. The vertices 
denote the associated breaking scales. At each of these breaking scales we need to apply the matching 



condition from Equation D.7 for the inverse gauge couplings. 



We find the following bound on the string scale, looking at the difference between the SU (3) 
gauge coupling and the U(1)y gauge coupling at Mz- 

a Y \M z ) - a^(M z ) = 89.88 < ^oT\M b ) + * log 

^M s > 10 14 GeV , (3.26) 
where in the last estimate we neglect the influence of the term involving a~ 1 (M s ). The 



details of this estimate can be found in Appendix DA . We remind the reader that the 
unification scale M s does not have to be the string scale but can be the string scale enhanced 
by the radius of the overall compactification M\jy ~ MgV 1 ^ 6 [34,35]. The above upper 
bound refers to the unification scale Mjjv an d hence the upper bound on the string scale is 
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lower. We can achieve for a string scale at the intermediate scale which depending on the 
structure of soft masses [21,22] can be consistent with TeV soft-masses. 

Using the running derived in Appendix [D], we find the following notable scenarios: 

1. Intermediate string scale: Breaking the copies of SU(A) and SU(2)r at the string 
scale, we find a scenario with gauge coupling unification, which saturates the lower 
bound on the string scale (cf. Equation 3.26| ) of 10 14 GeV. Note that for now the 
additional Higgs generations are allowed to survive until the weak scale. Bounds on 
flavour changing neutral currents generally forbid these but only to the PeV scale [68] , 
changing the nature of the running not significantly in this approximation. The 
running is illustrated in Figure 0. 
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Figure 6: Shown is the running for the inverse gauge couplings applicable in between the relevant 
scales. The green, blue and red circles at Mz denote the experimentally observed values for the 
gauge coupling constants. The orange circle at M a denotes the inverse gauge coupling at the string 
scale. 



2. GUT scale string scale: This is an example of a string scale that is close to the 
GUT scale at M s = 10 16 GeV. The breakdown to the MSSM gauge group is at the 
high scale. The running is illustrated in Figure |7[ 
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Figure 7: In similar fashion to the intermediate string scale scenario we show the running for the 
inverse gauge couplings. In the scenario on the left, the breaking down to the Standard Model gauge 
groups takes place at two stages, whereas the breaking in the scenario on the right takes place at 
the string scale at 10 16 and we observe the standard gauge coupling unification of the MSSM. 



3. Additional left-handed W bosons at the LHC: The breaking of the two SU(2)l 
gauge groups can be as low as the TeV scale from the perspective of gauge coupling 
unification, which offer interesting phenomenological opportunities in the LHC era. 
The running is visualised below in Figure ||. 



string 




x [GeV] 



Figure 8: Shown is the running for the inverse gauge couplings applicable in between the relevant 
scales. The green, blue and red circles at Mz denote the experimentally observed values for the 
gauge coupling constants. The orange circle at M s denotes the inverse gauge coupling at the string 
scale. 



4. Consistency of scales 

After discussing the model, we now comment on the consistency of scales associated to the 
breakdown to the Standard Model gauge symmetries and gauge coupling unification with 
the scales of moduli stabilisation and supersymmetry breaking. In particular we would like 
to compare it with the two supersymmetry breaking scenarios in the context of Large volume 
compactifications [15], depending on whether moduli redefinitions occur [21,22] or not [20]. 
Hierarchies are created by different suppressions in terms of the overall bulk volume in string 
units V. In both scenarios we control contributions to the scalar potential up to 1/V 3 . Hence 
we demand D- and F- flatness up to that order. We have to ensure that any contribution 
from the matter fields is sub-leading, and discuss in turn the requirements. The different 
scenarios for soft-masses essentially fix the overall bulk volume to provide a 'solution' to 
the hierarchy problem. In scenario 1 with soft masses of order m so f t ~ l/47r 2 Mp/V we find 
a volume V ~ 10 12 . In scenario 2 the smallest soft masses scale as m so f t ~ Mp/V 2 which 
then corresponds to a volume of order V ~ 10 7 . 
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These soft-masses are obtained when starting with the following [43,44] tree- level scal- 
ing in the matter Kahler potential 

K ™ = ^<h4, (4-1) 

where fa denotes any matter field. We satisfy the constraint from the contribution to the 
potential by requiring the fields to obtain vevs less than of the order 

M P 

<^ C > < ym ■ ^ 

In the largest case, the F-term of the matter field \F^\ = {K^F^F^ 1 / 2 scales as \F^\ ~ 
Mp/V 3 / 2 . 



In scenario 1, this limiting case of a suppression as in Equation if! is consistent with 
the smallness of the matter F-terms as required by the other soft-masses. Using a volume 
of order V ~ 10 12 , we hence have a maximal breaking scale of order 

(fa) ~ M b ~ 10 9 GeV. (4.3) 

In this setup, unification shall occur at M uv ~ Mp/V 1 / 3 ~ 10 14 GeV. At the superficial 
level with the exclusion of threshold effects and supersymmetry breaking effects in the 
running, we find the breaking at 10 9 GeV too small for gauge coupling unification. We note 
that the inclusion of warping effects or tuning Wq might evade this constraint. 

In the second scenario, the F-terms of the matter fields have to be suppressed up to 
1/V 2 and hence require the matter fields to obtain a vev of the order 

(fa) ~ M b ~ — ^ ~ 10 12 GeV, (4.4) 

where we assumed a volume of the order V ~ 10 6 . As shown in Figure |8| we find this limiting 
scale to be consistent with gauge coupling unification. 



5. Conclusions and Outlook 

We have studied the construction of realistic models on del Pezzo singularities and argued 
that in order to account for the hierarchies of masses, flavour structure and symmetry 
breaking patterns dPo,dP\ and (IP2 are not rich enough whereas dP^ satisfies all the re- 
quirements. We presented a brane model based on the third del Pezzo singularity dP% that 
gives rise to all Standard Model fermions with three pairs of Higgs-doublets at low-energies 
with realistic Yukawa structure for both quarks and leptons. The control over the flavour 
sector can arise through the control of the leading order Kahler potential for the matter 
fields, utilising the additional gauge structure of anomalous U(l) symmetries in dP^. These 
U(l) symmetries along with U(1)b-l play a crucial role in forbidding proton decay. To 
achieve the desirable flavour parameters (12 masses and 6 angles) we tune 14 parameters. 
We use a vev for the right-handed sneutrino to generate Majorana neutrino masses and 
checked that R-parity violation be sufficiently suppressed. Our model also predicts that the 
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angle #13 of the PMNS mixing matrix is non-zero but detailed values may be affected by 
loop effects. 

Depending on the scale of breakdown to the Standard Model gauge groups we found 
gauge coupling unification in the energy range from Mjjv ~ 10 14 GeV (which can corre- 
spond to a string scale at M s ~ 10 12 GeV) to the usual unification scale at Mtjv ~ 10 16 GeV, 
without the requirement of one GUT group at high energies. A priori, the lower bound on 
the string scale sounds very interesting in the context of standard soft-masses [20] being 
corrected by loop-effects to the order of the gravitino mass divided by loop suppression 
factors [21,22] since this still can lead to TeV scale soft-masses. Unfortunately, the break- 
ing scale of the additional gauge symmetries to the Standard Model is not unrelated to 
supersymmetry breaking since F-flatness is only guaranteed to the order the potential is 
controlled. This upper bound on the breaking scales renders it impossible to achieve gauge 
coupling unification, without the inclusion of threshold effects and supersymmetry break- 
ing effects, at an intermediate string scale unless tuning W$ or introducing warping effects. 
However, the upper bound on the breaking scale is found to be consistent with a string 
scale near the usual unification scale at 10 16 GeV and TeV soft-masses as in the scenario 
presented in [20] . A proper discussion of supersymmetry breaking terms is beyond the scope 
of this article. 

These detailed properties distinguish the model on (IP3 from previously presented mod- 
els as for example in [12] and [14]. In particular we should list the realistic Lepton Yukawa 
couplings, realising hypercharge from non-abelian factors hence guaranteeing it to be mass- 
less, the complete scenarios for flavour physics and unification presented as successes of the 
construction on (IP3. 

Depending on the breaking scales and the favourite model of unification, we find very 
interesting phenomenological signatures of our model, naively in the reach of LHC physics, 
such as additional U(l) symmetries, additional S'C/(2)-bosons and an interesting Higgs 
sector. 

However, at this stage of model building there still remain some open problems offering 
rich prospects for future work: 

• To achieve the desirable flavour and mass properties, we need to tune several vevs by 
hand. It remains open to determine the potential for these vevs and to ask whether 
these vevs actually can be achieved dynamically. 

• To explain via a concrete model, possibly in a concrete compact or semi-local setup, 
how the .D3 — D7 states get masses via D7 — D7 states such that they are not present 
at low energies. 

• To find an explicit realisation of the scenario with intermediate scale string scale and 
unification which required warping or tuning of the flux parameter Wo. 

• To obtain the correct scale for neutrino masses dynamically. 

• To find an embedding of this model on dP% in an honest-to-God string compactifica- 
tion in terms of a compact Calabi-Yau with dP% singularities. 
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The list of achievements and open questions illustrates that the bar on realistic string 
models is increasing with time which is encouraging. We hope to return to some of these 
questions in future work. 
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A. Consistency of Breaking 
A.l D-Flat 

After vevving D3 — D3 fields such as 3?6i to break down to the Standard Model, we must 
ensure that the supersymmetric D-terms are satisfied. The abelian D-term equations can 
always be satisfied by tuning the FI terms and they do not require further consideration. 
The non-abelian D-terms must be satisfied at each node of the quiver, and this requires 
that also D3 — D7 fields obtain vevs. To demonstrate this, consider the non-abelian D-term 
conditions at the SU(4) node, letting t a be the SU(4) generators. If we denote the incoming 
fields as X p and the outgoing fields as Y q , the following condition must be satisfied: 



£ X p t a X p = Y q t a Y Q . (A.l) 



Similar equations hold at all other nodes of the quiver. If no fields were vevved, Equation [AT 
would certainly be satisfied in the vacuum of the theory. It should also be clear that if an 
incoming field has a vev, then an outgoing field will also need a vev to satisfy the equation. 
As all fields in the quiver are bi-fundamentals, this vevved outgoing field will be the ingoing 
field for another node of the quiver whose D-term condition must be satisfied. Continuing 
to apply this reasoning around the quiver, we see that for the D-terms to be satisfied the 
vevved fields must form closed cycles in the quiver, and hence correspond to gauge invariant 
operators [69,70]. We can therefore satisfy all the D-term equations by vevving 37 string 
states associated with the vevved D3 — D3 state. 

One can easily check that the D-flat conditions for the VEV (ps3)i are satisfied but 
they are not satisfied for (^53)2 (all are satisfied apart from T 15 ). An additional VEV is 
needed which for example is given by a pair of D3 — D1 states, call them p 7a and p^. The 



same applies to (3>6i)- I n particular the vev in Equation |C.2| and |C,3| can be used to satisfy 



the non-abelian D-term equations, leading to conditions that can easily be satisfied 

= -H 2 + M 2 + |p| 2 , and = \cf>\ 2 - \4>\ 2 + \<p 7 \ 2 . (A.2) 
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A. 2 F-flat 

To satisfy D-flatness, we had to vev D3 — D7 states. Now one can check that the cubic 
coupling among the D3 — D3 and the D3 — D7 states is no longer F-flat in the global 
supersymmetric sense d%W 7^ 0. However we are working in an effective supergravity frame- 
work, implying that we should look at the associated F-term in supergravity meaning 
F l = K lJ DjW 7^ 0, where D^W = djW + WdiK denotes the usual covariant derivative in 
supergravity. In the effective supergravity we have only control over the (next-to)-leading 
order contribution to the potential and in this framework we can trust, respectively have 
to guarantee, F-flatness up to the order we can trust our effective supergravity theory. For 
concreteness we work in a large volume effective supergravity [15] where the expansion pa- 
rameter is the overall bulk volume and beyond the leading order contribution to the scalar 
potential, corrections in g s and a' ruin the approximation. The effective supergravity setup 
is given by 

K = -2 log (V + O + ^jlh: *i*J j , (A.3) 

W = W + Ae~ aT ° Wetter , (A.4) 

where V denotes the overall volume of the bulk geometry, £ the leading order a' corrections 
and Tj 4-cycle volumes, 3>j denotes all matter fields, Wo the flux parameter. For further 
details we refer the reader to the review literature on the subject [71,72]. In this supergravity 
setup we can show that depending on the size of the vev for the matter fields, we can achieve 
F-flatness. Assuming an overall scaling of ~ Mp/V ~ 1TI3/2, we find that all F-terms 
are suppressed at least at F % ~ I/V 7 / 3 . This is an additional suppression compared to the 
non- vanishing F-terms of the Kahler moduli by more than 1/V which for us is sufficient 
for F-flatness. Having a larger suppression for the vev of the field would render the F-term 
even smaller. 

A.3 Masses 

A vev for any bi-fundamental field can induce masses via the D-term couplings. Taking for 



example the diagonal vev for ^42 as in Equation |3.2| the induced masses are proportional to 
m tp ~ 9^- Out of the original eight real degrees of freedom in ^42, all of them apart from 
four, which are eaten as Goldstone bosons of the broken U(2), obtain a mass. The same 
can be checked explicitly in the SU (4) case. The D3 — D7 states do not obtain a mass via 
the U(2) D-term but obtain a mass via couplings with D7 — D7 couplings which are not 
part of the local construction. We conclude that after breaking all degrees of freedom in 
^42 become massive as desired. In a similar way it can be shown that all degrees of freedom 
involved in the breaking of the SU(4) and SU(2)r become massive. 

B. Non-anomalous U{\) symmetries 

As for example discussed in [73], the non-anomalous U(l) symmetries can be identified as 
follows: The anomaly of interest is U(l)i x SU(Nj) 2 . The anomalies are best summarised 
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in an n x n matrix whose entry (i, j) corresponds to the anomaly between U (l)j x SU (Nj) 
The entry T{j is given by 



m N 



fe=i 1=1 

where M 1 ^ 1 ) is a matrix that includes the £7(1) charges for all fields, M su ^ is the matrix 
with all charges for SU(N) gauge groups, $ indicates the multiplicity of each field and 
5(.) returns 1 only if that matrix element is non-zero. Generally one then finds that all 
U(l) symmetries will be anomalous for itself but that there are linear combinations of U (1) 
symmetries which will be anomaly-free. Here we are interested in linear combinations that 
involve only combinations of D3 gauge groups since in general the D7 gauge groups could 
have further charged matter which is not localised at the singularity, which influences the 
charges. Generally, one obtains one anomaly-free combination generated by 

Q anomaly— free — ^ ^ TV"' C^'^) 



For our Pati-Salam model on dP% it turns out by applying equation ( |B.1| ) that this is the 
only anomaly-free combination. This U(l) combination will however not turn out to be 
important and we assume that this cycle is non-trivial in a global embedding and have 
hence no massless U (1) factor from the abelian symmetries. 

C. Majorana masses via D-brane instantons 

We can wrap an E3 brane on either the cycle associated with m,2, or m\2 (cf. Figure^). 
We demand that they have the same Chan-Paton factors as the D3 branes so that they 
can be interpreted as gauge instantons. The arguments are similar to the above discussion 
of /U— term and we directly jump to the discussion of the charged zero modes. 

The cancellation of anomalies does not allow for a solution with all three cycles mi , 
or mi2 wrapped by E3 branes but only a combination of two of them, one example shown 
in Figure [|. However in this example we get a non- perturb at ive coupling for the remaining 
neutrino via the higher dimensional non-perturbative coupling 4a — 2Ra — 2Kb — E3 — 4a. 
It can be argued that such higher dimensional operators are present [74]. 

Furthermore to generate a quadratic coupling, such as the Majorana neutrino masses, 
we have to have two zero modes coupling to ur, which can be achieved by SP(2) instantons 
[50, 75]. In addition the unequal ranks (4 and 2) leading to unequal zero modes for the 
E3 — D3 states implies that it is necessary to have D7 — E3 zero modes to generate a non- 
vanishing contribution to the superpotential after integrating over the fermionic zero modes. 
We therefore require states between E3a, E3b and vti\\ fj , fj 2 where couplings including f/ 2 
involve higher dimensional couplings including ^53. Let us study an example which should 
be taken just showing the possibility of generating the coupling. The couplings present for 
the instanton zero modes are shown 

ViQfvl r)Mp73, ViQ§vf, ViViP53p73, *ji*37^\ vf^f, (C.l) 
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where i = 1,2. To give a contribution proportional to v R v R we ^ n( ^ constraints on the 
allowed rank for the D7 gauge groups. In our notation gauge group m\\ has to have rank 3 
and gauge group has to have rank 1. Other constructions involving D7 states coupling to 
Standard Model fields are dangerous since a vev would induce large unobserved masses for 
quarks. Furthermore we observe that the vev structure for D7 states $37 and p^ requires 
a breaking of SU{4) -> SU(3) x £7(1) and SU(2) R ->■ 17(1) : 



P73= * , (C2) 




$37 = , , (C3) 



where <p and p denote some vev which has to be in accordance with the D-flat condition. 
With this vev structure one can show that after integrating over the charged zero modes 



according to Equation 2.7 the only quadratic coupling induced in the effective action is 
the Majorana neutrino mass. We find that all three Majorana masses are distinct due to 
different 33 states appearing in the coupling ^53 and <E>61- Schematically we then have: 

^Majorana = A(j>, ^"^^l + ^f A (P, <\>Y~ a ^V 2 V 2 + ^A(p, ^e'^V^ , (C.4) 

where Tj denotes the chiral superfield associated to the 4— cycle volume the E3 brane is 
wrapping. Besides the dependence on 73 states in A, we expect A not to have any further 
suppression. 



However, for these masses to be present, the vev structure in Equation |C.2| and C.3 
requires that after the breaking to SU{3) x U(1)b-l an d U(V)2R these gauge groups are 
the diagonal product of the D3 and D7 gauge groups £7(3) or respectively U(l). For the 
U(l) symmetries to remain massless it is necessary that the abelian £7(1) factors in the D7 
gauge groups remain massless. This is not the case in our setup since the Kahler modulus 
is charged under these £7(1) symmetries and hence generates a Stiickelberg mass for these 
£7(1) symmetries. In this setup it is hence impossible to keep the standard hypercharge in 
models based on the Pati-Salam gauge groups massless. 

In principle, the hierarchy between the //—term (electroweak scale) and the Majorana 
mass (intermediate scale) can arise since they originate from distinct E3 branes wrapping 
different cycles. We shall note that all bi- linear operators in the right-handed quark fields 
are induced with the same suppression involving different D3 — D7 fields compared to the 
right handed neutrino bi-linear. We assume that the D3 — D7 fields do not obtain a vev 
and due to the large suppression of that coupling they are phenomenologically irrelevant. 

D. Details on gauge coupling unification 

In this appendix we give the detailed conventions and expressions used for the gauge 
coupling running and bother about details on the £7(1) normalisation factors. We take 
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the normalisation of the generators of the fundamental representation of U(N) to be 
Tr(T a T b ) = l/25 ab . In this convention we have for the SU(2) generators T a = where a a 
are the standard Pauli matrices. The diagonal generator in 577(4) T 15 leading to U(1)b-l 
is given by 

/:k \ 




-.15 



-h 
-h 



(D.l) 



We have the following Lagrangian for the gauge field strength and a field ip transforming 
in the fundamental of U(N) 



C = --F« V F»™ + ^ V (dp + igA«T a ) V> . 



(D.2) 



SU(N) a x SU(N) b -»• SU(N) d 

We start with the following Lagrangian for the bosonic components of the gauge fields and 
fields ^1,2 in the fundamental of SU(N) at b 

£ = \F%F»™ - \f%F^ + 4^ (d, + i 9a AlT a ) Y>i + 4^ (d, + igbB^) ^ 2 (D.3) 
Let us define the new gauge fields after breaking in terms of the old generators 



cos6> -sine \ I L°\ 
sm6 cos# J I M» ) ' 



(D.4) 



where denotes the generator of the diagonal gauge group and the orthogonal massive 
SU(N) generators. With this definition we see that the Lagrangian becomes 



£ = -/UK a - - A F* vM F^ b + 4^ (d, + iga [cos (9)L; - sin (9)M*] T a ) ^ 



% + "ft [sin {0)L% + cos (0)M£] T a ) V> 2 
\*%,L*T - \ F Ui F M h + rfV (d, + ig[Ll - tan (9)M*) T a ) ^ 
+4^ (dp + ig[L« + cot (6)M*\ T a ) ^ , 



(D.5) 



where we used the definition of the diagonal gauge coupling in terms of the old gauge 
coupling 

gd = g a cos 9, g d = g b sin 9 . (D.6) 

Independent on the breaking scale we find for the inverse gauge couplings the following 
boundary condition 
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SU(4) SU(3) x U(1) B -l 

We start with a Lagrangian for the bosonic gauge field and a fermionic field ip in the 
fundamental of SU{4) 



(D.8) 



After breakdown to SU(3) x U(1)b-l the surviving generators are the usual Gell-Mann 
generators and the diagonal generator T 15 (cf. Equation p.l| ). ift decomposes into q trans- 
forming as (3, 1/3) and e (1, —1). We hence can write after breaking 



1 



QB-h)q 



-(F 3 f - -(F 15 ) 2 + gty*^ + igA Al\ a + ig 4 A™T 15 )q 
+eV(^ + ig A A«\ a - i 94 Al\l -Qb-l U 



where we defined 



This directly leads to 



9B—L 



754 



a B-L ~ g a 4 



(D.9) 
(D.10) 
(D.ll) 



SU(2) R x U(1) B -l -> U(1) y 

Recall that ?7(l)y is generated by = 1/2Qb~l + This leads to the conventional 
charges with the left-handed quarks having charge 1/6. Starting with the following La- 
grangian 

--(F B . L f - \{F 2R f + (S M + 2ig B „ L A»-- L 



C = -^(F B -l) 2 - L A {F2 R f + 4^(9^ + igB-LA^QB-L^i + ^(fy + igmA^T )^ 



+4^{d tl + ig 2R Ai«T a )^ 



2Rrpa\ 



(D.12) 



In analogy to the discussion of the breaking of the diagonal gauge groups we then find gy 
in terms of the previous gauge couplings to be given by 



2g B -L sin ip = g y , g 2 R cos tp = g y 



This leads to 



ay 1 = a 2 \ + ^a B 1 _ L . 



(D.13) 
(D.14) 
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SU(2) L xU(l) Y 4U(l) em 

Electromagnetic U(l) then is generated by Q em = Qy + T§. In the same fashion as before 
in the case of breaking two gauge groups to the diagonal gauge group, one finds for the 
gauge coupling of electromagnetism e in terms of gy and g 2 i as follows 

e = gy cos 6w, e = g 2 Lsin#w. (D.15) 

This leads to the boundary conditions for the inverse gauge couplings at Mz 

(D.16) 



a em — a Y + a 2L 



Summary of gauge coupling running in dP^ 

This is a summary of the running of gauge couplings in the model on (IP3. For the original 
gauge groups it becomes 

= <*-\M s ) + ^ log , for M 4d < x < M s (D.17) 

^L,b(*) = *~\ M s) + %^ log (^f) , for M 2Ld <x<M s (D.18) 

a 2Ra,b( x ) = <*~\M S ) + log (^j , for M 2Rd <x<M s (D.19) 
We find the following (3— function coefficients 

/3 4a = (-3/2 + me, -7), (3 4b = (-9/2, -7) , (D.20) 
/3 2 La = (m 6 ,-l), /3 2Lfe = (4 + m 6 -m 9 ,0), (D.21) 
/?2*a = (5/2 + m 6 , 0), /5 2ra = (-1/2, -1) , (D.22) 

where the first number denotes the gauge coupling including D7 states in the running and 
the second one involves no D7 states in the running. After breaking to the diagonal gauge 
groups we have the following running for the gauge coupling running 

M 4d 



a 4d 1 (x) = ^ 1 (M 4d ) + ^log( 

= 2a'\M s ) + log (Jfy + ^ log (^) , for M 4 < x < M 4d (D.23) 

c£(x) = a-\M 2Ld ) + ^log (^) (D.24) 

= 2a- 1 (M s ) + ^L±^L log + ^ l og , for M z < * < M 2Ld 

a 2 l R (x) = a'\M 2Rd ) + ^ log (D.25) 

o 1 &Ra + &2Rb , ( M s \ (3 2Rd ( M 2Rd \ 

= 2a (Ms) + log — + — — log , for M z < x < M 2Rd 

2vr \M 2Rd J 2vr \ x J 
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Assuming all D7 states being integrated out, we find the following /3— function coefficients 

= (-6, -6) , (D.26) 
PiLd = (3, 1) , (D.27) 
fold = (3, 1) , (D.28) 

where the second number corresponds to the (3— function with only one Higgs generation 
as opposed to three. After breaking of the SU(A) — > SU(3) x U(V)b-L a t M4 we have the 
following running for the gauge couplings 



a B-L 



i L (x) = c£ L m + ^ log ffi = + log (^) , for My < X < M 4 



16 _! 4(/3 4a + M , / Ms \ , 4/3 4d , 

y a (Ms) + — 3^ — log Uw + ^ g 



M id \ + ^L lQg (M ± )(D 3Q) 



M 4 / 2vr \ x 



where we used the boundary condition found in Equation D.ll for the running of U(V)b-l- 
We find the following (3— function coefficients 

fa = -3 , Pb-l = 16 . (D.31) 

After breaking SU{2) R x U{1) B -L -> tf(l)y, where Qy = Tf + Qb-l/2, we have the 
following running for the inverse gauge coupling for U(1)y 

a Y l = a^(Mr) + ^a£_ L (M Y ) + 7^ log (^) , for M z < x < My 

= H a -i (M- ) + ^ + ^ log (J^-) + ^ log f **™\ + log f "A 

3 1 s; 2^ 8 VAf2W 2tt B UW 3tt 8 V^W 

+ 3^ log {-m) + ^r log {M^) + ^ log J • (D - 32) 

For a MSSM matter content, the f3— function should be the standard one of 11 whereas 
with three Higgs generations it is 13 

(3 Y = (11, 13). (D.33) 

Note that in the case when we break to the MSSM gauge group at the string scale (inte- 
grating out the additional Higgs generations) the running reduces precisely to the running 
of the gauge couplings in the MSSM with the following boundary conditions 



tg' 1 = a 2 1 = 2a i 1 and a Y l = — a i 1 . (D.34) 
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D.l Bound on the string scale 

To achieve the experimentally observed difference between the SU(3) and U(l)y gauge 
coupling at the weak scale, we find a lower bound on the possible string scale in (IP3. This 
limit is ought to change when taking threshold and supersymmetry breaking effects into 
account. However, using the gauge coupling running from the previous sections for now we 
find 

-1 -1 W ^ B 2 Ra+Jhm. ( M s \ p 2Rd ( M 2Rd \ 

+ y±±M log (Mi) + *i log f^i) + *t± log (M±) + * i og (Ml) 

3tt \M m J 3tt & V M 4/ 8vr \M Y J 2tt \M z J 

- V a {Ms) + ~^r- log [wj + ^ log ) + log {m- z ) ) 

= 4 -a~HM s ) + + log f + ^ log - ^L±M log f Mi) 

log (M") + ^ log W + Al log f ^ _ A log (D.35) 

We now distinguish the following cases, involving first running with D3 — D7 brane states 
and in the second without. In both cases the difference is larger when we include all three 
Higgs generations. 

Estimate with D3 — D7 states 

To maximise the difference we find the following two cases depending on the values for m% 

Case A: p 2Ra + 8 2Rb > 8 2Rd , (2 + m 6 > 3) M 2Rd = My (D.36) 
Case B: p 2Ra + B 2Rb < 8 2Rd , (2 + m 6 < 3) M 2Rd = M s (D.37) 
Bia + fab >fad, (-5 + m 6 > -6) =► M 4d = M s (D.38) 



In the first limit Equation D.35 can be rewritten as 

ay 1 - a 3 - 1 = ^a-\M s ) + log (M s ) (^^) + log (M 4 ) 



+ ]o K ( Uv )(^p) -log(A^)^ 



<^(M s) + log(^)0) (D.39) 

where in the last line we maximised the difference by taking the breaking scales M4 and 
My to their maximal value M s . Similarly in the second limit we observe 

ay 1 -Q3 1 = ~ar l {M s ) +log(M s ) (A) +log(M 4 ) (A 



+ log (My) f — \ -log (M Z ) ( ~ 
<^- 1(M .) + 1 o g (^)(i) (D.40) 
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where the upper limit appears again in the limit My and M 4 approaching the string scale 



Estimate without D3 — D7 states 

In this case we observe 



P2Ra + p2Rb>p2Rd, ("K 3) M 2Rd = M s (D.41) 
/?4a + fab < fad, ("14 < -6) M 4d = M 4 (D.42) 



In this limit Equation D.35 can be rewritten as 



ay 1 -"* 1 = V 1 (M s )+log(M s ) (^^) + log (M 4 ) 
+ log (My) -log(M z ) 

In consistency with the breaking pattern, we find in both cases the same lower bound 
on the string scale. 

Estimate of a' 1 (M s ) 

The above estimate contains a term proportional to a~ 1 (M s ) and having in mind that we 
are able to tune a~ l (M s ) it is interesting to estimate its size given the experimental input 
of the gauge couplings at Mz- One straight forward constraint can be obtained from the 



gauge coupling of SU(2)l as given in Equation D.25 



2S.57 = a -i(M z) = + log (-^) + % !o g (^M) 

where in the last step our choice of beta functions is such that it would give the largest 
negative contribution due to the running, hence allowing the maximal value for a (M s ). 
We hence find the following upper bound on a~ 1 (M s ) 

a' 1 (M s ) < 17.36. (D.45) 

This upper bound is achieved in the limit M 2 id Mz and no D7 states in the running. 
Furthermore this estimate now enables us to give a lower bound on the possible string 
respectively breaking scales by combining it with the above estimate as discussed in the 
main text. 
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E. Pati-Salam without D7 branes 



The model discussed in the main part of the text requires a large D7 sector due to anomaly 
cancellation. One might wonder whether a different choice of D3 gauge groups can evade 
this bound. Figure ||| shows a C/(4) 4 x U(2) 2 gauge theory on dPs that satisfies the anomaly 
cancellation condition with unequal gauge group ranks. This is an extension of the Pati- 
Salam model in the main section in the left-right U(2) factors. For a realistic model the 
breakdown to the Standard Model is crucial and it turns out, using the techniques available 
to date, that it requires the re-introduction of D7 branes. We hence at this stage do not 
pursue this option any further. 




Figure 9: A Pati-Salam model on dPz which does not require D7 branes from anomaly cancellation. 
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